Superfluid Bose gas in two dimensions 
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We investigate Bose-Einstein condensation for ultracold bosonic atoms in two-dimensional sys- 
tems. The functional renormalization group for the average action allows us to follow the effective 
interactions from molecular scales (microphysics) to the characteristic extension of the probe I 
(macrophysics) . In two dimensions the scale dependence of the dimensionless interaction strength 
A is logarithmic. Furthermore, for large I the frequency dependence of the inverse propagator be- 
comes quadratic. We find an upper bound for A, and for large A substantial deviations from the 
Bogoliubov results for the condensate depletion, the dispersion relation and the sound velocity. The 
melting of the condensate above the critical temperature Tc is associated to a phase transition of 
the Kosterlitz-Thouless type. The critical temperature in units of the density, Tc/n, vanishes for 
/ — > 00 logarithmically. 
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I. INTRODUCTION 

Bosc-Einstein condensation and superfluidity for cold 
nonrclativistic atoms can be experimentally investigated 
in systems of various dimensions Two dimensional 
systems can be achieved by building asymmetric traps, 
resulting in different characteristic sizes for one "trans- 
verse extension" It and two "longitudinal extensions" I 
of the atom cloud [2]. For I ^ It the system behaves 
effectively two-dimensional for all modes with momenta 
(f' < l^^. From the two-dimensional point of view. It sets 
the length scale for microphysics - it may be as small as a 
characteristic molecular scale. On the other hand, the ef- 
fective size of the probe I sets the scale for macrophysics, 
in particular for the thermodynamic observables. 

Two-dimensional superfluidity shows particular fea- 
tures. In the vacuum, the interaction strength A is di- 
mensionless such that the scale dependence of A is log- 
arithmic [3']. The Bogoliubov theory with a fixed small 
A predicts at zero temperature a divergence of the occu- 
pation numbers for small q — \q\, n{q) ^ nc S^'^\q) Q. 
In the infinite volume limit, a nonvanishing condensate 
nc = Po is allowed only for T = 0, while it must van- 
ish for r > due to the Mermin- Wagner theorem Q. 
On the other hand, one expects a critical temperature Tc 
where the superfluid density po jumps by a finite amount 
according to the behavior for a Kosterlitz-Thouless phase 
transition [6i|. We will see that Tc/n (with n the atom- 
density) vanishes in the infinite volume limit ^ — > oo. 
Experimentally, however, a Bose-Einstein condensate can 
be observed for temperatures below a nonvanishing criti- 
cal temperature Tc - at first sight in contradiction to the 
theoretical predictions for the infinite volume limit. 

A resolution of these puzzles is related to the simple 
observation that for all practical purposes the macro- 
scopic size I remains finite. Typically, there will be a 
dependence of the characteristic dimensionless quantities 
as po/n, Tc/n or A on the scale I. This dependence is 
only logarithmic. While X{n ~ T ^ 0,1 ^ oo) = 0, 
{po/n){T / 0,/ oo) 0, {Tc/n){l ^ 0) = 0, in ac- 



cordance with general theorems, even a large finite I still 
leads to nonzero values of these quantities, as observed 
in experiment. 

An appropriate tool for studying the scale dependence 
of various quantities is the functional renormalization for 
the average action Q. Within a two-dimensional con- 
text, this starts with a given microphysical or classical 
action at some ultraviolet momentum scale Auv ~ Z^^. 
Then a scale dependent effective action F^ is defined by 
integrating out successively the modes with momenta \q\ 
larger than some infrared cutoff scale k. The flow of Ffc 
with k obeys an exact renormalization group equation, 
which we solve approximately by truncation of the most 
general form of F^. When k reaches a scale fcph ~ , all 
fluctuations are included since no larger wavelength are 
present in a finite size system. The experimentally rele- 
vant quantities and the dependence on / can be obtained 
from Ffcpj,. 

This paper investigates the functional renormalization 
group for a nonrelativistic Bose gas in two spatial dimen- 
sions. A priori, this method is not restricted to a pertur- 
bative treatment for a small coupling A and can therefore 
be used to investigate the effects of fluctuations beyond 
the Bogoliubov approximation. We assume that the in- 
teraction can be approximated by a contact term. This 
should be appropriate especially for dilute gases where 
the details of the interaction are not relevant. More ex- 
plicitly our microscopic action reads 



(1) 



We work with the functional integral within the Matsub- 
ara formalism. The "Euclidean time" r parameterizes a 
torus with circumference T~^, and the bosonic flelds 4> 
are complex, depending on the two-dimensional coordi- 
nates X and T, X = {x,t). In addition to the tempera- 
ture T and the chemical potential fi the action involves 
the dimensionless coupling A as a parameter. We use 
h ^ ks = 1 and units of energy where the mass of the 
atoms is set to 1/2, 2M = 1. 
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All macroscopic properties following from the micro- 
scopic model ([1]) can be obtained from the full effec- 
tive action T[(j)\. For example, the thermodynamic prop- 
erties can be inferred from the grand canonical parti- 
tion function Z or the corresponding grand potential 
= —TlnZ. The latter is related to the effective ac- 
tion by 



(2) 



where the expectation value 0oq is obtained as a solution 
of the field equation 



In 



0. 



(3) 



Since the effective action is the generating functional 
of the one particle irreducible (IPI) correlation functions 
we can also derive dynamical properties from F [</>]. For 
a homogeneous field (j>cq the full propagator G{p) is ob- 
tained from 



S(j)*{pi) 5(j),(j)2) 



= {G-%ipi)d{pi~p2). 



(4) 



Here we employ a basis with two real fields, i = 1,2, 

0=75(01+^02), 0*(p) =0.(-p). 

The calculation of the functional r[(j)] is not an easy 
task - it includes the summation of infinitely many di- 
agrams in perturbation theory. Our method determines 
F[0] by using an exact flow equation for the scale depen- 
dent effective action Ffc[0]- the average action. Due to 
the presence of an infrared cutoff at the scale k only the 
quantum and thermal fluctuations with momenta ^ k'^ 
are included. In the infinite volume limit F^ interpolates 
between the microscopic action and the full effective ac- 
tion 
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For a system with finite size I we are interested in F^^j^, 
^ph — l^^ ■ If statistical quantities for finite size systems 
depend only weakly on I, they can be evaluated from 
Ffepj^ in the same way as their thermodynamic infinite 
volume limit follows from F. Details of the geometry etc. 
essentially concern the appropriate factor between /cph 
and . 

The average action obeys the exact flow equation 0, Q 



9..F,[0] = iTr(F(^)[0] 



Rkr^dkRk 



(6) 



Despite its simple one loop structure it includes all orders 
in perturbation theory as well as nonperturbative effects. 
Nevertheless, in praxis this equation is not solvable ana- 
lytically since it is a differential equation for a functional. 
Approximate solutions can be found by truncating the 



space of possible functionals Tk[4>\ to a manageable size. 
Such a truncation does not have to rely on an expansion 
in a small parameter, as it it the case for perturbation 
theory. Because of this, the method is nonperturbative 
and also applicable for large interaction strength. 

In this work we use a very simple truncation which 
includes terms with up to two derivatives 

+2^(m-Mo)0*(5.- A)0 + C7(p,Ai)l. (7) 



Here the effective potential U contains no derivatives and 
is a function of p = <t>*<t>- At the end we will evaluate all 
quantities for a flxed chemical potential /io ~ the term 
~ (/i — /io) will be employed for an easy evaluation of 
9Ffc/9/i, as needed for the density. In addition to the cou- 
plings parametrizing U (see below) our truncation con- 
tains three further /c-dependent ( "running" ) couplings A, 
S and V . The particular form of our truncation is mo- 
tivated by a more systematic derivative expansion and 
an analysis of symmetry constraints (Ward identities) in 
Ref. [9] . (This reference may also be consulted for details 
of the formalism and the flow equations.) A truncation 
similar to Eq. (17]) has been investigated at zero temper- 
ature in Refs. [lol . ITl| . A simpler truncation was used 
in Ref. [T^. In terms of renormalized flelds cj) = A^l'^cj), 
p = Ap, renormalized kinetic coefficients S — V — ^ 
and effective potential U{p,p) — U{p,fi), Eq. (O reads 

Tk - J^S^r{Sdr-A-Vd^)cp 

+2V{^l-^lo)(|)*{^r-A)cf, + u{p,^i)}. (8) 



For the effective potential, we use an expansion around 
the fc-dependent location of the minimum po{k), and the 
fc-independent value of the chemical potential pQ that 
corresponds to the physical particle number density n. 
We determine poik) and po by the requirements 

{dpU){po{k),po) -0 foraUfc 
-{df,U){po, Pa) =n at fc = fcph. (9) 

More explicitly, we take a truncation for U{p,p) of the 
form 

U{p,p) = U{po,po) -nk{p- po) 

+ [m^ + a{p - po)] [p - po) 

+ l[X + fi{p-p,)]{p-p,)\ (10) 

In the symmetric phase we have po = 0, while in the 
phase with spontaneous symmetry breaking, we have 
rri^ = 0. The atom density n = —dU /dp corresponds 
to rik in the limit k — > fcph. 



In total, we have at fixed chemical potential /i = /io 
and fixed temperature T four running renormalized cou- 
plings po{k), A(fc), S{k), and V{k). (In the symmetric 
phase po is replaced by m^.) In addition, we need the 
anomalous dimension 77 — — fc9fclnA. To determine the 
flow of the density nk, we also need the dependence of 
the effective potential on (/i — /xq) and we therefore in- 
clude the couplings a{k), f3{k). The precise procedure 
how the flow equation for the average action is projected 
to flow equations for the running couplings is explained 
in Ref. 0. The flow equation for the effective potential 
Uk {p, p) and the kinetic coefficients S and V as well as 
the anomalous dimension rj are given there for arbitrary 
dimension d. In this paper we apply these flow equations 
for d = 2. 



II. SCALE DEPENDENCE OF INTERACTION 
STRENGTH 
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FIG. 1: (Color online) Flow of the interaction strength A(fc) 
at zero temperature and density for different initial values 
Aa = 100, Aa = 10, Aa = 4, Aa = 2, Aa = 1, and Aa = 0.4 
(from top to bottom). 



It is well known that in two dimensions the scattering 
properties cannot be determined by a scattering length 
as it is the case in three dimensions. In experiments 
where a tightly conflning harmonic potential restricts the 
dynamics of a bose gas to two dimensions, the interac- 
tion strength has a logarithmic energy dependence in the 
two-dimensional regime [12]. For low energies and in the 
limit of vanishing momentum the scattering amplitude 
vanishes. In our formalism this is reflected by the loga- 
rithmic running of the interaction strength A(fc) in the 
vacuum, where both temperature and density vanish. 
In general, the flow equations in vacuum describe the 
physics of few particles like for example the scattering 
properties or binding energies. Following the calculation 
of Ref. , we flnd the flow equations for the interaction 
strength (t = ln(fc/A)) 

A2 (fc2 - m2) 

Since in vacuum the propagator is not renormalized, 
dtS = dtV = dtA = dtm^ = 0, we set S* = 1 and F = 
on the right-hand side of Eq. (|lip . The vacuum corre- 
sponds to — [9] and we obtain the flow equation 



(12) 



The vacuum flow is purely driven by quantum fluctua- 
tions. It will be modifled by the thermal fluctuations for 
T 7^ and for nonzero density n. 
The solution of Eq. (fT2l) 



A(fc) = 



1 



TI + 4^1n(A/fc) 



(13) 



goes to zero logarithmically for fc ^ 0, X{k = 0) = 0. 
In contrast to the three-dimensional system, the flow of 
the interaction strength A does not stop in two dimen- 
sions. To relate the microscopic parameter Aa to exper- 
iments exploring the scattering properties, we have to 



choose a momentum scale Qoxp, where experiments are 
performed. To a good approximation the relevant in- 
teraction strength can be computed from Eq. by 
setting k = Qoxp- If not specifled otherwise, we will use a 
renormalized coupling A = A(fcph). 

For our calculation we also have to use a microscopic 
scale A below which our approximation of an effectively 
two dimensional theory with pointlike interaction be- 
comes valid. Our two-dimensional computation only in- 
cludes the effect of fluctuations with momenta smaller 
then A. In experiments A~^ is usually given either by the 
range of the van der Waals interaction or by the length 
scale of the potential that conflnes the system to two 
dimensions. We choose in the following 



A = 10, fcph = 10 



(14) 



At this stage the momentum or length units are arbitrary, 
but we will later often choose the density to be n = 1, 
so that we measure length effectively in units of the in- 
terparticle spacing n^^^^. For typical experiments with 
ultracold bosonic alkali atoms one has n"^/^ « lO^^cm. 

The flow of A(fc) for different initial values Aa is shown 
in Fig. [1] Following the flow from A to kph yields the 
dependence of A = A(/cph) on Aa as displayed in Fig. [2l 
It follows from Eq. (|13p that for positive initial values 
Aa the interaction strength A is bounded by 



A < 



47r 



ln(A/A:ph) 



1.82. 



(15) 



The last equation holds for our choice of A and fcph. We 
emphasize that our bound holds only if the interactions 
are approximately pointlike for all momenta below A. 
Close to a Feshbach resonance this may not be true and 
our formalism would need to be extended by considering 
nonlocal interactions or introducing an additional field 
for the exchanged two-atom state in the "closed channel" . 
In other words, close to a Feshbach resonance the effective 
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FIG. 2: (Color online) Interaction strength A at the macro- 
scopic scale feph = 10~^ in dependence on the microscopic 
interaction strength Aa at A = 10 (solid). The upper bound 
-^max = in(A/fcp,,) is also shown (dashed). 



A, 5, V 











1.0 






/ 

/ 

/ 

^ 0.8 

/ 

/ 

/ 






! 0.6 

/ 






/ 

/ 

/ 0.4 

/ 

/ 






" ' -A . 0.2 

















-10 -8 -6 -4-2 2 



FIG. 3: (Color online) Flow of the kinetic coefficients A 
(solid), S (dashed), and V (dashed-dotted) at zero tempera- 
ture r = 0, density n = 1, and vacuum interaction strength 
A = 1. 



cutoff A for the validity of a two-dimensional model with 
pointlike interactions may be substantially lower, and the 
upper bound on A correspondingly higher. 

III. FLOW EQUATIONS AT NONZERO 
DENSITY 

In this section we investigate the many body problem 
for a nonvanishing density n at zero temperature T = 0. 
A crucial new ingredient as compared to the vacuum is 
the nonzero superfluid density 

"-s = Po(fcph)- (16) 

For interacting bosons at zero temperature the density n 
and the superfluid density ng are equal. In contrast, the 
condensate density is given by the unrenormalized order 
parameter 

nc = Paikph) = ^"^(fcph)po(fcph)- (17) 

Due to the repulsive interaction, nc may be smaller than 
n, the difference n — nc being the condensate depletion. 
(In the limit A — > we have n = ns = nc-) To obtain 
a nonzero density at temperature T = we have to go 
to positive chemical potential /i > 0. At the microscopic 
scale k — A the minimum of the effective potential U is 
then at po,A = m/^ > 0. 

The superfluid density po is connected to a nonvanish- 
ing "renormalized order parameter" 00, with po = 4iq4io- 
It is responsible for an effective spontaneous breaking of 
the U(l)-symmetry. Indeed, the expectation value (j)o 
points out a direction in the complex plane so that the 
global U(l)-symmetry of phase rotations is broken by the 
ground state of the system. Goldstones theorem implies 
the presence of a gapless Goldstone mode, and the as- 
sociated linear dispersion relation uj ^ \q\ accounts for 
superfluidity. The Goldstone physics is best described 



by using a real basis in field space by decomposing the 
complex field = 0o + -^(0i +«</'2)- Without loss of gen- 
erality we can choose the expectation value (j>o to be real. 
The real fields (pi and 02 describe then the radial and 
Goldstone mode, respectively. For p = po the inverse 
propagator reads in our truncation 

_ T ff + VpI + U' + 2pU'\ ^Spo \ 
^ ~ [ Spo, pF + Vql + U')- 

(18) 

Here p is the momentum of the collective excitation, and 
for r = the frequency obeys oj — —ipa- In the regime 
with spontaneous symmetry breaking, po{k) ^ 0, the 
propagator for p = po has U' = 0, 2pU" ~ 2Xpo ^ 0, 
giving rise to the linear dispersion relation characteristic 
for superfluidity. This strongly modifies the flow equa- 
tions as compared to the vacuum flow equations once 

<C 2Apo- For n ^ the flow is typically in the regime 
with po{k) ^ 0. In practice, we have to adapt the ini- 
tial value po,A such that the flow ends at a given density 
/Oo(^ph) = n. For fcph <C n^/^ one finds that po(^ph) de- 
pends only very little on fcph. As mentioned above, we 
will often choose the density to be unity such that ef- 
fectively all length scales are measured in units of the 
interparticle distance n~^/^. 

In contrast to the vacuum with T = pQ — 0, the flow 
of the propagator is nontrivial in the phase with po > 
and spontaneous U(l) symmetry breaking. In Fig. Owe 
show the flow of the kinetic coefficients A, V, S for a 
renormalized or macroscopic interaction strength A = 1. 
The wavefunction renormalization A increases only a lit- 
tle at scales where k « n^/^ and saturates then to a value 
A> I. As will be explained below, we can directly infer 
the condensate depletion from the value of A at macro- 
scopic scales. The coefficient of the linear r-derivative 
S goes to zero for fc — > 0. The frequency dependence is 
then governed by the quadratic r-derivative with coeffi- 
cient V, which is generated by the flow and saturates to 
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FIG. 4; (Color online) Flow of the interaction strength A(fc) 
at zero temperature T = 0, density n = 1, for different initial 
values Aa. The dotted lines are the corresponding graphs in 
the vacuum n = 0. The vertical line labels our choice of 
fcph. The lower plot shows X{k)/k for the same parameters, 
demonstrating that A(fc) ~ k for small k. 



a finite value for fc ^ 0. 

The flow of the interaction strength X{k) for different 
values of A = A(A:ph) is shown in Fig. ID While the 
decrease with the scale k is only logarithmic in vacuum, 
it becomes now linear X{k) ~ k for k ^ n^/^. It is 
interesting that the ratio X{k)/k reaches larger values for 
smaller values of Aa- 

As k is lowered from A to fcph, the renormalized order 
parameter or the superfluid density po increases first and 
then saturates to po = = 1- This is expected since 
the superfluid density equals the total density at zero 
temperature. In contrast, the bare order parameter pq 
Pq/A flows to a smaller value po < po. As argued in 
the bare order parameter is just the condensate density, 
such that 



n-nc = Pa- Pa= Po(l - -j) 



(19) 



is the condensate depletion. Its dependence on the inter- 
action strength A is shown in Fig. [51 For small interaction 
strength A the condensate depletion follows roughly the 
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FIG. 5: (Color online) Condensate depletion (n — n^jn as a 
function of the vacuum interaction strength A. The dashed 
line is the Bogoliubov result in — n^jn = ^ for reference. 



Bogoliubov form 



X_ 
Stt' 



(20) 



However, we find small deviations due to the running of 
A which is absent in Bogoliubov theory. For large inter- 
action strength A ~ 1 the deviation from the Bogoliubov 
result is quite substantial, since the running of A with the 
scale k is more important. 

We also investigate the dispersion relation at zero tem- 
perature. The dispersion relation w(p) follows from the 
condition 



detG"i(w(p),p) = 



(21) 



where is the inverse propagator after analytic con- 
tinuation to real time po As was shown in Q the 
generation of the kinetic coefficient V by the flow leads 
to the emergence of a second branch of solutions of Eq. 
(I^Tj) . In our truncation the dispersion relation for the 
two branches (p) and w_ (p) are 



1 

^(p' + Vo) + ^ 



l/2\ 1/2 



-— p^(p^ + 2Apo) 



(22) 



In the limit — > 0, S* — > 1 we find that the lower branch 
approaches the Bogoliubov result uj^ ^/p^p^ + 2Xpo) 
while the upper branch diverges uj+ — > oo and thus dis- 
appears from the spectrum. The lower branch is domi- 
nated by phase changes (Goldstone mode), while the up- 
per branch reflects waves in the size of po (radial mode). 

In principle, the coupling constants on the right-hand 
side of Eq. (|22p also depend on the momentum p = |p| . 
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FIG. 6: (Color online) Lower branch of the dispersion relation 
cj_(p) at temperature T — and for the vacuum interaction 
strength A = 1 (solid curve), A — 0.5 (upper dashed curve), 
and A = 0.1 (lower dashed curve). The units are set by the 
density n = 1. We also show the Bogoliubov result for A = 1 
(upper dotted curve) and A — 0.5 (lower dotted curve). For 
A = 0.1 the Bogoliubov result is identical to our result within 
the plot resolution. 
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FIG. 7: (Color online) Dispersion relation uj-{p), a;+(p) at 
temperature T — and for vacuum interaction strength A = 1 
(solid), A = 0.5 (long dashed), and A = 0.1 (short dashed). 
The units are set by the density n — 1. 



Since an external momentum provides an infrared cutoff 
of order k ~ {p] we can approximate the |p|-dependence 
by using on the right-hand side of Eq. the k- 

dependent couphngs with the identification k — \p\. Our 
result for the lower branch of the dispersion relation is 
shown in Fig. [6l W e also plot the Bogoliubov result 
uj = y/p2^p2 _|_ 2Xpo) for comparison. For small A our 
result is in agreement with the Bogoliubov result, while 
we find substantial deviations for large A. Both branches 
u)^ and CO- are shown in Fig. [7] on a logarithmic scale. 
Since we start with V = a.t the microscopic scale A we 
find uj+ip) oo for \p\ A. 

The sound velocity cs can be extracted from the dis- 
persion relation. More precisely, we compute the mi- 
croscopic sound velocity for the lower branch lo~{'p) as 
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FIG. 8: (Color online) Dimensionless sound velocity cs /n^/^ 
as a function of the vacuum interaction strength (solid). We 
also show the Bogoliubov result cs = \/2Apo for reference 
(dashed) . 



cs ~ ^ &t p = 0. In our truncation we find @ 
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(23) 



Our result for C5 at T = is shown in Fig. [5] as a 
function of the interaction strength A. For a large range 
of small A we find good agreement with the Bogoliubov 
result c|. = 2Apo. However, for large A or result for cg 
exceeds the Bogoliubov result by a factor up to 2. 



IV. NONZERO TEMPERATURE 

At nonzero temperature and for infinite volume, long 
range order is forbidden in two spatial dimensions by 
the Mermin- Wagner theorem. Because of that, no 
proper Bose-Einstein condensation is possible in a two- 
dimensional homogeneous Bose gas at nonvanishing tem- 
perature. However, even if the order parameter vanishes 
in the thermodynamic limit of infinite volume, one still 
finds a nonzero superfluid density for low enough tem- 
perature. The superfluid density can be considered as 
the square of a renormalized order parameter po = \4'q\'^ 
and the particular features of the low-temperature phase 
can be well understood by the physics of the Goldstone 
boson for a p hase with effective spontaneous symmetry 
breaking [l3|. The renormalized order parameter <j)Q is 
related to the expectation value of the bosonic field (^0 
and therefore to the condensate density po = by a 
wave function renormalization, defined by the behavior 
of the bare propagator G at zero frequency for vanishing 
momentum 

00=^1/2^0, Po^Apo, G-i(p^O)=A/. (24) 

While the renormalized order parameter pQ{k) remains 
nonzero for k — > if T < Tc, the condensate density 
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Po = Po/-^ vanishes since A diverges with the anoma- 
lous dimension, A ~ . After restoring dimensions the 
relation 



Po 



Po 

hm „ -, ■ 



(25) 



is the Josephson relation |14l |. 

The strict distinction between a zero Bose-Einstein 
condensate po = and a nonzero superfluid density 
Po > for nonzero temperature < T < Tc is valid only 
in the infinite volume limit of a homogeneous system. For 
a finite size of the system, as atoms in a trap, the running 
of A{k) is effectively stopped at some scale kph- There 
are simply no collective modes with wavelength larger 
than the size of the system, whose fluctuations would be 
responsible for a further increase of A. With a finite 
both Po and po are nonzero for T < Tc, and the distinc- 
tion between a Bose-Einstein condensate and superfluid- 
ity is no longer relevant in practice. For large systems 
^(^ph) can be large, however, such that the condensate 
density can be suppressed substantially as compared to 
the superfluid density. In any case, there is only one 
critical temperature Tc, defined by po{T < Tc) > 0. 

The flow equations permit a straightforward compu- 
tation of po{T) for arbitrary T, once the interaction 
strength of the system has been fixed at zero temperature 
and density. We have extracted the critical temperature 
as a function of A = A(/cph) for different values of fcph- 
The behavior for small A, 



An 



n ln(C/A) 



(26) 



is compatible with the free theory where Tc vanishes for 
kph — * and with the perturbative analysis in Ref. ^15'|. 
We find that the value of C depends on the choice of 
fcph. For fcph — 10^^ we find ( — 100, while fcph = lO^'' 
corresponds to C = 225 and fcph = 10"'^ to C = 424. 
In Fig. [9] we show or result for Tc/n as a function of 
A for these choices. We also plot the curve in Eq. 
with the Monte-Carlo result ^ = 380 from Ref. [161] . 
We find that Tc vanishes for fcph ^ in the interacting 
theory as well. This is due to the increase of ( and, for a 
fixed microscopical interaction, to the decrease of A(fcph). 
Since the vanishing of Tc/n is only logarithmic in fcph, a 
phase transition can be observed in practice. We find 
agreement with Monte-Carlo results |l6l| for small A if 
fcph /A « 10""^. The dependence of Tc/ n on the size of 
the system fcpj^^ remains to be established for the Monte- 
Carlo computations. 

The critical behavior of the system is governed by a 
Kosterlitz-Thouless phase transition. Usually this is de- 
scribed by considering the thermodynamics of vortices. 
In Refs. [13, [3 it was shown that functional renor- 
malization group can account for this "nonperturbative" 
physics without explicitly taking vortices into account. 
The correlation length in the low-temperature phase is 
infinite. In our picture, this arises due to the presence of 
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FIG. 9: (Color online) Critical temperature Tc/n as a func- 
tion of the interaction strength A. We choose here fcph — 10~^ 
(circles), fcph — lO"* (boxes) and fcph = 10~® (diamonds). 
For the last case the bound on the scattering length is 



A < 



0.78. We also show the curve — — 



ln(A/fcph) ..^ "..ov. o..„,. ^ - i„(C/A) 

(dashed) with the Monte-Carlo result = 380 [l^ for refer- 



a Goldstone mode if po > 0- The system is superfluid for 
T < Tc- The power like decay of the correlation function 
at zero frequency 



G{p) - if) 



(27) 



is directly related to the running of A. As long as fc^ 3> 
the bare propagator obeys approximately 



G 



1 



A{k)p 



Aik) 



(28) 



Once fc^ ^ p^, the effective infrared cutoff is given 
by instead of fc^, and therefore A{k) gets replaced 
by A{y/^), turning Eq. ^ into Eq. For 
large po the anomalous dimension depends on po and 
T, 77 = T/(47rpo). 

Another characteristic feature of the Kosterlitz- 
Thouless phase transition is a jump in the superfluid den- 
sity at the critical temperature. However, a true discon- 
tinuity arises only in the thermodynamic limit of infinite 
volume (fcph 0), while for finite systems (fcph > 0) 
the transition is smoothened. In order to see the jump, 
as well as essential scaling for T approaching Tc from 
above, our truncation is insufficient. These features be- 
come visible only in extended truncations that we will 
briefly describe next. 

,2 

For very small scales -^^ ^ 1, the contribution of Mat- 
subara modes with frequency qo = 2T:Tn, n ^ 0, is sup- 
pressed since nonzero Matsubara frequencies act as an in- 
frared cutoff. In this limit a dimensionally reduced theory 
becomes valid. The long distance physics is dominated 
by classical two-dimensional statistics, and the time di- 
mension parametrized by r no longer plays a role. 

The flow equations simplify considerably if only the 
zero Matsubara frequency is included, and one can use 
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FIG. 10: (Color online) Flow of the density n (solid), the 
superfluid density po (dashed), and the condensate density 
Po (dotted) for chemical potential p = 1, vacuum interaction 
strength A = 0.5 and temperatures T = Q (top), T — 2.4 
(middle) and T = 2.8 (bottom). The vertical line marks our 
choice of kph- We recall n = po for T = such that the upper 
dashed and solid lines coincide. 

more involved truncations. Such an improved truncation 
is indeed needed to account for the jump in the super- 
fluid density. In Ref. [l^ the next to leading order in a 
systematic derivative expansion was investigated. It was 
found that for fc <C T the flow equation for po can be well 
approximated by 

dtpo = 2.54r-i/2(o.248T-po)^^^fi'(0.248r-po). (29) 

We switch from the flow equation in our more simple 
truncation to the improved flow equation (|29p for scales 
k with fc^/T < 10"'^. We keep all other flow equations 
unchanged. A similar procedure was also used in Ref. 

In Fig. [TO] we show the flow of the density n, the 
superfluid density po and the condensate density po for 
different temperatures. In Fig. [TT]we plot our result for 
the superfluid fraction of the density as a function of the 
temperature for different scales fcph- One can see that 
with the improved truncation the jump in the superfluid 
density is indeed found in the limit kph 0. Fig. [T2l 
shows the condensate fraction po/n and the superfluid 
density fraction po/n as a function of T/n. We observe 
the substantial /cph dependence of the condensate frac- 
tion, as well as an effective jump at Tc for small fcph. We 
recall that the infinite volume limit fcph = amounts to 
Po = for T > 0. 

The Kosterlitz-Thouless description is only valid if the 
zero Matsubara frequency mode (n = 0) dominates. For 
a given nonzero T this is always the case if the the char- 
acteristic length scale goes to infinity. In the infinite 
volume limit the characteristic length scale is given by 
the correlation length ^. The description in terms of a 
classical two dimensional system with U(l) symmetry is 
the key ingredient of the Kosterlitz-Thouless description 
and holds for ^^T 1. In the infinite volume limit this 



Po/n 




FIG. 11: (Color online) Superfluid fraction of the density po/n 
as a function of the dimensionless temperature T/n for inter- 
action strength A = 0.5 at different macroscopic scales fcph = 1 
(upper curve), fcph = 10"°'^, fcph = 10~\ fcph = 10~^-^, 
fcph = lO"'^, fcph = lO"'^'^ (bottom curve). We plot the result 
obtained with the improved truncation for small scales (solid) 
as well as the result obtained with our more simple trunca- 
tion (dotted). (The truncations differ only for the three lowest 
lines.) 



po/n, po/n 




FIG. 12: (Color online) Condensate fraction of the density 
Po/n as a function of the dimensionless temperature T/n for 
interaction strength A = 0.5 at macroscopic scale fcph = 10"'^ 
(solid curve) and fcph = 10~* (dashed curve). For compari- 
son, we also plot the superfluid density po/n at fcph = 10^^ 
(dotted). These results are obtained with the improved trun- 
cation. 



always holds for T < Tc or near the phase transition, 
where ^ diverges or is very large. For a finite size system 
the relevant length scale becomes fcpj^^ if this is smaller 
than ^. Thus the Kosterlitz-Thouless picture holds only 
for T > kl^. 

For very small temperatures T < fcpj^ one expects a 
crossover to the characteristic behavior near a quantum 
critical phase transition, governed by the quantum criti- 
cal fixed point. The crossover between the different char- 
acteristic behaviors for T > fcpj^ and T < k^^ can be 
observed in several quantities. As an example we may 
take Fig. [TO] and compare the flow of po and po for low 



9 



T (close to the T = curve) or large T (other curves). 



CONCLUSIONS 



is dominated by phase fluctuations and shows the 
characteristic superfluid behavior w p for p — > 0. 
The upper branch corresponds to excitations of the 
local superfluid density po- 



In conclusion, we have demonstrated that functional 
renormalization for the average action yields a unified 
picture for ultracold nonrelativistic bosons in arbitrary 
dimension d. We have employed in this paper the same 
simple truncation for d = 2 as in Ref. [gj for d = 3. While 
for d = 3 the perturbative Bogoliubov treatment is reli- 
able for a wide range of parameters, we find for d — 2 
substantial deviations if the coupling A is large, or for 
the regime close to the phase transition. The Kosterlitz- 
Thouless type of the transition is nonperturbative in A 
and cannot be seen within the Bogoliubov approxima- 
tion. 

In particular, we have obtained the following main re- 
sults for interacting Bose gases in two dimensions: 

1. An upper bound on the interaction strength A in 
vacuum arises if the characteristic momentum of 
the scattering particle p = |p| is sufficiently be- 
low the ultraviolet scale A below which a two- 
dimensional description with pointlike interactions 
becomes valid. 



\{p)< 



ln(A/p) ■ 



(30) 



The dependence of A(p) on p is logarithmic, 
d\/d\Dp = AV(47r). 

2. The condensate depletion at T = is reasonably 
well approximated by the Bogoliubov result if A is 
small, (n — nc)/n = A/(87r). Logarithmic correc- 
tions arise due to the running of A. For large A 
the repulsive interaction suppresses the condensate 
and the condensate depletion can exceed the Bo- 
goliubov result by a large factor. 

3. For the dispersion relation at T = we observe 
a substantial deviation from the Bogoliubov result 
to = p\/p'^ -\- 2Xn for large A. We also find a second 
branch of "radial excitations". The lower branch 



4. The sound velocity at T = exceeds the Bogoli- 
ubov result for large A by up to a factor 2, while 
the perturbative result is valid for small A. 

5. The system shows a second order phase transition. 
The critical temperature Tc depends logarithmi- 
cally on A, Tc/n = 47r/ln(C/A). The coefficient C 
depends logarithmically on the size of the system 
I = fc-j^i. For /cph « lO^^A we find C ~ 400 and 
agreement with Monte-Carlo simulations. The l- 
dependence of Tc/n still needs to be established by 
numerical simulations. 

6. The phase transition is found to be of the 
Kosterlitz-Thouless type, with a jump of the super- 
fluid density Apo ~ Tc/4:. For a flnite size system 
this jump is smoothened. Below Tc, the correlation 
length is inflnite. The correlation function shows a 
powerlike decay, with an anomalous dimension 77, 
G{p) ^ p^"^^^ . The anomalous dimension depends 
on the superfluid density po ^md on the tempera- 
ture rj = T/(47rpo)- 

7. We have computed the superfluid density po a-nd 
the condensate density nc = po as a function of 
T. The condensate density vanishes in the infinite 
volume limit / ^ cx), if T > 0. For finite I a nonva- 
nishing condensate can be observed. 

In this paper we have employed a rather simple trun- 
cation which is sufficient to show all the characteristic 
features of two-dimensional cold and dilute nonrelativis- 
tic bosons. One exception is the jump in the superfluid 
density which needs an improved truncation in the range 
of very small fc^/T. Extensions of the truncation are 
straightforward and should improve the quantitative ac- 
curacy for large A. With increasing experimental preci- 
sion the verification of a quantitatively precise theoretical 
computation will be a challenge. 
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